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The combination of electronic correlation and spin-orbit coupling is thought to precipitate a va¬ 
riety of highly unusual electronic phases in solids, including topological and quantum spin liquid 
states. We report a Raman scattering study that provides evidence for unconventional excitations 
in a-RuCla, a spin-orbit coupled Mott insulator on the honeycomb lattice. In particular, our mea¬ 
surements reveal unusual magnetic scattering, typified by a broad continuum. The temperature 
dependence of this continuum is evident over a large scale compared to the magnetic ordering tem¬ 
perature, suggestive of frustrated magnetic interactions. This is confirmed through an analysis of 
the phonon linewidths, which show a related anomaly due to spin-phonon coupling. These obser¬ 
vations are in line with theoretical expectations for the Heisenberg-Kitaev model and suggest that 
a-RuCla may be close to a quantum spin liquid ground state. 


The quantum spin liquid, where quantum fluctuations 
obstruct long-range order at low temperatures, is one of 
the most elusive and intriguing states of matter[l, 2]. 
These states have been suggested to play a crucial role 
in high-temperature superconductivity[3] and could lead 
to topologically protected excitations useful for quantum 
computation[4]. The elementary excitations of these sys¬ 
tems are fractionalized, reflecting the peculiar nature of 
the quantum spin liquid (QSL) ground state, and can be 
identihed in scattering experiments where they manifest 
themselves as broad continua, in contrast to the sharp 
magnon modes characteristic of ordered magnets [5-8]. 
While QSLs have been extensively studied theoretically, 
the experimental picture is limited due to the small num¬ 
ber of candidate material systems[l]. 

In this context, it was recently suggested[9, 10] that 
strongly spin-orbit coupled Mott insulators on the hon¬ 
eycomb lattice might manifest a two-dimensional QSL 
ground state. In such a material, strong spin-orbit 
coupling can lead to the formation of jeff = 1/2 
pseudospins [9-11], local moments of mixed spin and 
orbital character, whose bond-dependent exchange in¬ 
teractions map onto the Kitaev Hamiltonian: Hx = 
-JKj2-y-Hnks^i^]- Here Sj refers to the 7 compo¬ 
nent of the spin at the lattice site and the sums are 
performed over the three different types of links of the 
honeycomb lattice ^ = x,y, z. For a given site, the ex¬ 
change interactions along all three links cannot be si¬ 
multaneously satisfied. These bond-dependent interac¬ 
tions are therefore intrinsically frustrated, in a manner 
distinct from geometric frustration, and lead to an ex¬ 
otic, spin-disordered ground state and fractionalized Ma- 


jorana fermion excitations[12]. Experimental efforts in 
this direction have so far focused on the honeycomb lat¬ 
tice iridates[13-15], although these compounds are struc¬ 
turally complex [16, 17] and the applicability of Kitaev 
physics has also been questioned [17-19]. 

To experimentally explore spin liquid behavior driven 
by spin-orbit coupling, we have investigated the ele¬ 
mentary excitations of a-RuCla using polarized Raman 
scattering. a-RuCla is a recently identified spin-orbit 
assisted Mott insulator [20] which crystallizes in a lay¬ 
ered structure with planes of edge-sharing RuCle octa- 
hedra arranged in a honeycomb lattice[21] [Fig. 1 (a)] 
and is therefore an excellent candidate for realizing the 
Heisenberg-Kitaev model in the solid state[9, 10]. The 
magnetic properties of a-RuCla have been subject to a 
number of previous studies[21-23]. Magnetic suscepti¬ 
bility and specific heat measurements suggest a pair of 
magnetic phase transitions near 8 K and 14 K[24, 25]. 
Intriguingly, a recent single crystal neutron diffraction 
experiment suggests a zigzag or stripy-type magnetic or¬ 
dering with a small ordered moment[25], as expected 
for a Heisenberg-Kitaev magnet close to a QSL ground 
state[26, 27]. However, no spectroscopic information ex¬ 
ists regarding the spin excitations of a-RuCls and the 
possibility of spin fractionalization driven by spin-orbit 
coupling, a defining feature of the Kitaev QSL state, 
therefore remains an open question. Raman spectroscopy 
is an excellent tool for studying this possibility, as it 
can simultaneously probe spin, lattice, and electronic 
excitations[28, 29]. 

Our measurements reveal anomalous magnetic light 
scattering in a-RuCla, typihed by a broad continuum 
below 100 K that persists in the magnetically ordered 
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FIG. 1. Structure and polarized Raman response of a-RuCls. 

(a) Lattice structure of a single plane of a-RnCR. The Rn 
atoms are shown in blue while the Cl are indicated in red. 

(b) 5 K Raman intensity, (c) Low-energy detail at 5 K. The 
shaded blue region is a guide to the eye and indicates the 
continunm contribution, (d) 295 K Raman intensity. The 
Raman spectra evince a series of narrow phonon modes in 
both channels, a broad continuum extending up to 20 meV in 
Eg, and a quasielastic scattering component in Aig. The con¬ 
tinuum emerges at low temperatures, while the quasielastic 
scattering is only visible at high temperatures. 


state(s) below 14 K. This behavior is not easily explained 
by conventional two-magnon scattering or structural dis¬ 
order, but is consistent with theoretical predictions for 
the Kitaev spin liquid [8]. These results suggest that a- 
RuCls hosts unusual magnetic excitations and may be 
proximate to a QSL ground state. 

The a-RuCla crystals used in this study were grown 
by vacuum sublimation of prereacted RuCla powder [20]. 
The Raman spectra were measured in the quasi- 
backscattering geometry in both collinear (XX) and 
crossed (XY) polarizations, with light polarized in the 
basal (cleavage) plane. Light from a 532 nm laser was fo¬ 
cused down to a 10 /rm spot and the power at the sample 
is estimated to be 500 fiW. Two holographic notch filters 
(Ondax) were used to reject light from the fundamental, 
leading to a lower cutoff of 2 meV in XY and 3.5 meV 
in XX. The resolution of our spectrometer is estimated 
to be 0.6 meV. No anti-Stokes scattering was observed 
at 5 K, suggesting negligible laser-induced heating. This 
was confirmed by halving the power and doubling the in¬ 
tegration time, which produced no noticeable change in 
the spectra recorded at 5 K. 

The polarized Raman intensity of a-RuCla is shown 
in Fig. 1 for 5 and 295 K. The 5 K spectrum in Fig. 1 
(b) consists of a series of sharp phonon modes superim¬ 
posed on a weak continuum that extends up to roughly 
20 meV. The continuum is shown in detail in Fig. 1 (c). 
Strong and sharp phonon modes are resolved at 14, 20, 


34, 37, and 38 meV, with weaker modes detected at 27 
and 42 meV. Due to the broad line shape and low en¬ 
ergy, we attribute the continuum to magnetic scattering. 
In comparison, an electronic mechanism for the contin¬ 
uum is unlikely as a-RuCls is electrically insulating with 
a transport gap of 0.3 eV [30]. As we describe later, the 
magnetic assignment is further justified by the temper¬ 
ature dependence of the continuum scattering, which is 
inconsistent with a phonon mechanism. 

The phonon spectra can be understood in terms 
of an isolated layer of D^d symmetry, as anticipated 
from the layered crystal structure(see the Supplemental 
Material[31]. The observed modes are sharp, suggestive 
of a well-ordered crystal lattice, and show no evidence for 
structural transitions down to 5 K. We note that the two 
lowest energy phonons of Eg symmetry near 14 and 20 
meV show Fano line shapes at low temperatures [32]. This 
effect arises generically when a narrow resonance couples 
to a continuum. Furthermore, the Fano line shape of 
the 20 meV phonon demonstrates that the continuum 
extends up to ~20—25 meV. 

Having identified the excitations present, we focus on 
the temperature dependence of the magnetic scattering. 
For conventional (e.g. phonon) scattering, the Raman 
intensity should show a monotonic decrease given by the 
thermal Bose factor, although this is not the case in a- 
RuCls. In Fig. 2 (a), we plot the Eg Raman intensity 
at several temperatures. The continuum loses intensity 
before reaching a minimum near 100 K. Below 100 K, 
the continuum gains intensity without a change in en¬ 
ergy scale or shape, in contrast to conventional scatter¬ 
ing. This temperature dependence can be appreciated by 
considering the spectral weight SW = l{uj)duj, where 
the cutoffs uJi and U!h are 2.5 and 12.5 meV, shown in 
Fig. 2 (b). SW decreases before increasing below 100 K, 
at odds with the trend expected for one- or two-particle 
scattering due to the thermal Bose factors {[n(a;) -I- 1] 
or [n(a;) -|- l]^}(See the Supplemental Material[31]). The 
change in the magnetic scattering near 100 K is consis¬ 
tent with the suggested onset of in-plane spin correla¬ 
tions, first noted by Kobayashi et al.[22] in relation to 
magnetization measurements. Importantly, the persis¬ 
tence of magnetic scattering far above Tjv is a signature 
of frustrated, low dimensional magnetism. 

A comparison of the 100 K and 5 K curves in Fig. 2 (a) 
reveals that the continuum shape does not vary strongly, 
even below the ordering temperatures of 8 and 14 K. 
Typically, the broad continuua observed in the param¬ 
agnetic state of an antiferromagnet will evolve into well- 
defined peaks due to one or two magnon excitations[33]. 
In contrast, the scattering in a-RuCla displays no obvious 
changes upon ordering. We suspect that lower energies 
or temperatures than are accessible in our experiment 
may be needed to detect the magnons associated with 
the ordered magnetic state. 

The observed magnetic continuum is unusual and its 
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FIG. 2. Magnetic scattering in RuCls. (a) Raman intensity in Eg and A 2 g (XY). (b) Eg and A 2 g (XY) low-energy spectral 
weight (SW) vs. temperature, (c) Raman intensity in Aig (XX-XY). (d) A\g (XX-XY) low-energy spectral weight (SW) vs. 
temperature. As shown in (c), the Aig quasielastic scattering (QES) decreases in intensity as temperature is reduced down to 
100 K. In contrast, the intensity of the Eg continuum decreases as temperature is reduced to 100 K before increasing. This 
is more clearly resolved in the spectral weights shown in b) and d). The behavior of the Eg continuum is at odds with the 
temperature dependence of the thermal factors expected for one [n(a;) -|- 1] and two [(n(cu) -1- 1)^] particle scattering and signals 
a change in the magnetic excitations near 100 K. 


origin is not immediately clear. X-ray absorption and 
ab initio studies have firmly established the relevance of 
jeff = 1/2 physics[ 20 , 26], while static probes of the mag¬ 
netism have hinted at the existence of the Kitaev inter¬ 
action in q;-RuC 13 [ 24 , 25]. Scattering from fractionalized 
spin excitations is therefore an appealing interpretation 
for the magnetic continuum. Indeed, a recent theoreti¬ 
cal work identified a broad Eg continuum, extending up 
to 3 Jr- (where Jk is the Kitaev interaction as defined 
earlier), as the Raman signature of the Kitaev spin liq¬ 
uid [ 8 ]. Equating ZJk with the experimental upper cut¬ 
off of continuum scattering of 20—25 meV yields Jk 
8 meV. On a more general, empirical level, the broad 
line shape and temperature dependence of the contin¬ 
uum in a-RuCla are reminiscent of the behavior found in 
the putative spin liquid materials herbertsmithite[ 6 , 34] 
and Cs 2 CuCl 4 [ 7 ][ 35 ]. In these compounds, the contin¬ 
uum scattering gains intensity as spin liquid correlations 
develop. CS 2 CUCI 4 also orders at low temperatures. Sim¬ 
ilar to a-RuCla, the continuum scattering in CS 2 CUCI 4 
displays a minimal change upon entering the ordered 
state, with well-defined magnons only evident at very 
low energies and involving a small portion of the to¬ 
tal spectral weight [7]. Finally, a recent Raman study of 
(Nai_a;Lia;) 2 lr 03 also reported a broad continuum that 
was interpreted in terms of Kitaev QSL physics[36]. 

A second parallel between a-RuCl 3 and herbert- 
smithite is the presence of quasielastic scattering 
(QES) [ 6 ]. As can be seen in Fig. 1 (d), a QES com¬ 
ponent emerges in a-RuCls at high temperatures in the 
Aig channel. Such scattering is observed in a variety of 
low-dimensional magnetic systems [28] and is usually as¬ 
signed to fluctuations in the magnetic energy density. In 
Fig. 2 (c), we plot the Aig intensity - 4j/ 

for select temperatures. This QES intensity decreases 
down to 100 K, signifying a gradual evolution of the low- 


energy spin dynamics, at which point it becomes difficult 
to resolve in our data. This trend is borne out in the 
Aig spectral weight, integrated from 4 to 17 meV, shown 
in Fig. 2 (d)[37]. This could be due to a reduction in 
the intensity or scattering rate of the quasielastic fluc¬ 
tuations, but we cannot discriminate between these two 
possibilities with our experimental cutoff. Similar to the 
continuum scattering, the QES displays a strong temper¬ 
ature dependence far above the ordering temperature, as 
expected for a low-dimensional, frustrated magnet [28]. 

In systems with moderate spin-phonon coupling, the 
lattice dynamics provide an additional probe of the mag¬ 
netic degrees of freedom. Indeed, the Fano line shape of 
the Eg phonons is a signature of spin-phonon coupling 
in of-RuCls. The spectral region depicted in Fig. 3 (a) 
is well fitted at 5 K by a function consisting of the sum 
of two Fano line shapes and a constant background term 
(See Supplemental Material[31]). In Figs. 3 (b) and (c) 
we show the linewidth (T) and mode frequency (wo) for 
the 20 meV phonon. As shown in Fig. 3 (c), the linewidth 
r(T) displays a nonmonotonic temperature dependence. 
It decreases almost linearly with reduced temperature 
from 300 K down to 140 K. The width then begins to 
grow before plateauing at low temperature. 

Phonon self-energies are typically determined by lat¬ 
tice anharmonicity, leading to a monotonic temperature 
dependence [4]. However, r(T) reveals a nonmonotonic 
behavior. To demonstrate this, we have included in Figs. 
3 (b) and (c) plots of the behavior expected from an¬ 
harmonicity (See Supplemental Material[31]). r(T) and 
the fit diverge near 140 K, implying an additional relax¬ 
ation mechanism. This is seen in Fig. 3 (d), where we 
plot AE = r — Tank, which reveals an onset near 140 K, 
similar to the 100 K scale associated with anomalies in 
the continuum scattering and magnetic susceptibility[ 22 ]. 
We conclude that there are two contributions to the low- 
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FIG. 3. Spin-phonon coupling in a-RuCla. a) Data and Fano 
fit to the low-energy scattering at 5 K. b) Energy ujo of the 
20 meV phonon, c) Linewidth F of the 20 meV phonon, d) 
Magnetic contribution AF to the phonon linewidth. The low- 
energy data are well-described by the Fano form and indicate 
a coupling between the magnetic continuum and the lattice. 
This coupling is also apparent in the temperature dependence 
of r, which shows an anomaly near 140 K. 


lying Eg phonon self-energies: an anharmonic term and 
a component due to spin-phonon scattering which onsets 
near 140 K. This results in an increase in T and to the 
Fano line shape of the Eg phonons. The phonon proper¬ 
ties shown in Fig. 3 do not show obvious changes in the 
magnetically ordered phases below 14 K, suggesting that 
any effects due to spin-phonon coupling upon ordering 
are subtle. 

Overall, our experimental results appear consistent 
with the notion that a-RuCla hosts unusual magnetic 
excitations, possibly driven by proximity to a Kitaev 
QSL state. At this point, however, we consider other 
possible interpretations of the continuum, namely heav¬ 
ily phonon-damped two-magnon excitations and stacking 
disorder. Spin-phonon coupling is evident in the data and 
suggests an alternative explanation for the magnetic con¬ 
tinuum: conventional magnons that are heavily damped 
by magnon-phonon scattering. Indeed, for sym¬ 
metry, two-magnon (2M) scattering is also expected in 
the Eg channel[39]. Additionally, magnon-phonon cou¬ 
pling is known to cause unusually broad 2M features in 
some cases, for instance in NiPS3[39] and the insulating 
cuprates[40]. However, we regard the scenario of heavily 
damped 2M scattering as unlikely in a-RuCla for several 
reasons. First, the line shape of the magnetic scattering 
does not change appreciably when going from 100 to 5 K, 
in contrast to previous reports of phonon-damped 2M 
scattering[39, 40]. For dominant magnon-phonon scat¬ 
tering, the decay rate of a magnon of energy u) should 


be roughly proportional to the thermal occupation fac¬ 
tor n(a;)-|-l of a phonon of the same energy[40]. Tak¬ 
ing the magnon energy to be 10 meV (the center of the 
continuum), this would imply a 20% change in lifetime 
at low temperature. For a 5 meV magnon, this change 
would be 40%. Experimentally, however, no significant 
change is observed in the continuum scattering over this 
temperature range. Similarly, the characteristic energy 
of 2M scattering typically blueshifts considerably as tem¬ 
perature is lowered and evolves into a well-defined peak 
below Tff [39], in contrast to our observations in a-RuCla 
[Fig. 2 (b)]. We therefore conclude that spin-phonon 
coupling does not appreciably damp the magnetic exci¬ 
tations in a-RuCla. This is a result of the relatively low 
energy scale of the magnetic excitations compared to the 
phonons. In compounds where the magnons are heavily 
phonon damped, the characteristic energies of the mag¬ 
netic excitations are high compared to the phonon en¬ 
ergies, and so a large number of phonon branches are 
able to participate in the decay process[39]. In contrast, 
the magnetic excitations in o-RuCla are far less energetic 
than the majority of phonon branches present, meaning 
the phase space for spin-phonon scattering is significantly 
reduced. 

A further issue is the possible influence of stacking 
faults. Of-RuCla is a layered material and so such de¬ 
fects are likely present. Indeed, a short out of plane (c 
axis) magnetic correlation length was reported in a re¬ 
cent neutron diffraction study and attributed to stack¬ 
ing faults [25]. This stacking disorder in fact suggests 
another interpretation of the low-energy continuum re¬ 
ported here. In Raman scattering, stacking faults can 
lead to a relaxation of the k = 0 selection rule. In heavily 
disordered SiC, for example, this leads to both a broad¬ 
ening of existing Raman features and the appearance of 
new ones, including broad low-energy scattering corre¬ 
sponding to defect-activated acoustic phonon modes[41]. 
It is not clear, however, how static structural disorder 
could produce the unusual temperature dependences and 
QES present in our data. In the same vein, one might 
expect a broadening of the magnetic excitations due to 
stacking faults on the order of the interplane exchange 
interactions. Since the electronic structure of a-RuCla is 
highly two-dimensional[20], we expect that the interplane 
exchange interactions, and thus any broadening due to 
stacking faults, to be small compared to the overall en¬ 
ergy scale of the magnetic continuum, which is set by 
intraplane interactions. 

To conclude, we have studied the elementary excita¬ 
tions of a-RuCla using Raman scattering. Our measure¬ 
ments reveal a continuum of scattering that is not read¬ 
ily explained by conventional 2M scattering or structural 
disorder. Rather, the broad continuum reported here 
appears consistent with theoretical expectations for the 
Kitaev spin liquid. However, we stress that a categorical 
assignment is not currently possible and will likely require 
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inelastic neutron scattering. We observe changes in the 
magnetic dynamics (as evinced by the magnetic contin¬ 
uum and QES, as well as the 20 meV phonon self-energy) 
near 100—140 K, which we interpret as the development 
of in-plane spin correlations. Overall, our study suggests 
that a-RuCla may be close to a Kitaev spin liquid state, 
consistent with recent neutron diffraction and theoretical 
studies[25, 26]. 
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SUPPLEMENTAL MATERIAL 


I. Phonon Mode Assignment 


/(w) 


Is jq + e^) 

— 1 1 + 


( 1 ) 


Structurally, a-RuCla consists of a stack of weakly 
van der Waals bound layers. Accordingly, the number 
and symmetry of the phonons modes can be understood 
in terms of an isolated RuCla layer of symmetry D^d, 
rather than the full space group. In this case, group the¬ 
ory predicts a total of six (4ifg and 2Aig) Raman-active 
modes[l]. The AEg modes are expected to be visible in 
both polarization channels while the 2Aig modes vanish 
in the crossed geometry. Indeed, the modes at 38 and 
42 meV are suppressed in the crossed (XY) geometry 
and so can be assigned Aig symmetry. The remaining 
five modes as well as the continuum are present in both 
polarization channels and therefore have Eg symmetry, 
meaning RuCla exhibits one more Eg phonon mode than 
expected from group theory for a single layer. We 
therefore suspect that the mode at 27 meV is either de¬ 
fect activated or due to interlayer interactions, as it is the 
weakest peak in the Eg channel. Given that the spectra 
show only weak deviations from the ideal D^d symmetry 
of a single layer, we conclude that interlayer lattice in¬ 
teractions are indeed weak in this compound and so do 
not significantly affect the phonon dynamics. Although 
all observed phonon modes broaden in going from 5 to 
295 K, the number of modes does not change, meaning 
the D^d symmetry of the individual RuCla layers is un¬ 
perturbed. 


Here Is is the integrated intensity, q is the Eano asym¬ 
metry parameter, and e = 2(w — Wo)/r is a reduced 
frequency determined by the phonon frequency ujo and 
width r. The factor of — I in the denominator is a 
normalization factor chosen such that the integrated in¬ 
tensity is independent of q. 

As detailed in the main text, the low-frequency region 
is well-described by a Eano parameterization. A fit to 
a sum of three Lorentzians (one for the continuum and 
two for the phonons) was also tried, but yielded a con¬ 
sistently worse fit (e.g. the reduced was a factor of 
three larger than that of the Eano parameterization at 
low temperatures). 

In determining the temperature dependence of the 20 
meV mode parameters lOq and E, we limited our fit to 
the scattered intensity in the energy range within 2 meV 
of the phonon peak in order to avoid the complicating 
effects of the thermal Bose factor. The fact that the 
phonon frequency uio does not evince any deviations from 
the anharmonic form is in fact consistent with a con¬ 
tinuum density of states that does not vary strongly in 
the vicinity of uJq- This further confirms the presence 
of a broad magnetic continuum rather than well-defined 
magnon modes that are strongly damped by phonons. 
This leads to a constant (in energy) T and a correspond¬ 
ingly small frequency renormalization [2]. 


II. Thermal Bose Factor and Spectral Weight 

To estimate the temperature dependence of the scat¬ 
tering from purely thermal effects we adopt the following 
procedure. First, we divide the 295 K spectrum by a fac¬ 
tor of n{oj, 295 K) + l. We then multiply this quantity by 
n{uj,T) -I- I to generate the spectra anticipated at other 
temperatures. Finally, these spectra are integrated as 
described in the text to obtain the spectral weight (SW). 
For two-particle scattering, a factor of (n(w,T) -|- I)^ is 
used instead. 


III. Eano Resonance Fits 

The Eano effect can occur when a well-defined excita¬ 
tion, in this case an optical phonon, overlaps and also 
interacts with broad continuum of typically electronic 
or magnetic origin. These two excitations can interfere, 
leading to an asymmetric (Eano) resonance. The scatter¬ 
ing profile of this resonance as a function of frequency w 
is[2]: 


IV. Phonon Anharmonicity 

Phonon self-energies are typically determined by an¬ 
harmonic decay. Decay into a pair of acoustic modes of 
opposite momenta yields the following expressions for the 
phonon line width and frequency[3]: 

r(r) = r„ + A[i + 2n(w„/2)], ( 2 ) 

u;{T)=uJo-B[l + 2n{u;o/2)], (3) 

where tUo is the bare phonon frequency. To is the zero- 
temperature line width, A and B are constants, and 
n(^) is the Bose factor for a phonon of energy Wo/2. 
Importantly, equations 2 and 3 are monotonic functions 
of temperature, varying quasi-linearly for temperatures 
greater than cOo while changing slowly at temperatures 
less than oJq- In general, other combinations of acoustic 
and optical phonons that satisfy energy and momentum 
conservation can also contribute[4]. However, these pro¬ 
cesses lead to a similar monotonic temperature depen¬ 
dence that is almost constant at low temperatures and 
quasi-linear at high-temperatures. As discussed in the 
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main text, the temperature dependent line width of the 
20 meV Eg phonon cannot be described by this expres¬ 
sion. In the fits shown in the text, To and A were chosen 
from a fit to the high temperature linear regime between 
160 K and 300 K, while uJq was Hxed to the 5 K value 
derived from the Fano fits. 
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